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Abstract 

Analyses of QSO absorption lines are showing that HI content has evolved over the redshift 
range z=5 to z= 0. The 21cm line measurements of the z=0 HI content avoid several biases 
inherent in the absorption line technique, such as the influence of evolving dust content in 
the absorbers, and will produce a reliable measure to anchor theories of galaxy evolution. 
Examples of important questions to be addressed by local HI surveys are: (1) is there a 
significant population of gas-rich galaxies or intergalactic clouds that are missing from the 
census of optically selected galaxies? (2) is there an adequate reservoir of neutral gas to 
substantially prolong star formation at its present rate? and (3) are there massive objects 
of such low HI column density that they can have escaped detection in the “unbiased” HI 
surveys that have been conducted so far? 
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1 Evolution of the HI Content 

A measurement of the total HI content contained in nearby galaxies and HI clouds is an important 
constraint in the bigger picture of galaxy evolution on cosmological time scales. Although the 
neutral gas in large galaxies at present is often considered to be a minor component that is used 
as tracer for kinematics or as a dwindling source of fuel for star formation, there is now strong 
evidence from the studies of QSO absorption lines that the present HI content of the universe is but 
a fraction of what it was at redshift z ~ 3. This finding comes from the statistics of the “damped 
Lyman-ct” class of absorption line that is identified with dynamically cold layers akin to the disks 
of familiar nearby spiral galaxies (Wolfe et al 1986, Lanzetta et al 1995, Storrie-Lombardi et al 
1996). Figure 1 summarizes the presently available measurements of the neutral gas content as a 
function of z (Storrie-Lombardi et al 1996), along side a plot of the incidence of CIV absorption 
lines, which provide an indication of the cross section presented by clouds of ionized, metal-rich gas 
(Steidel 1990, Bahcall et al 1993), and the comoving density of luminous optically selected QSOs 
(Schmidt et al 1994, Hewett et al 1993). Taken literally, the evidence points to an epoch around 
z ~ 3 when the neutral gas mass density, Q g , hit a maximum, at roughly the same time that 
the CIV cross section for strong absorption lines (H / rest (A1548) > 0.15 A) began a sharp increase. 
This is also a time when luminous QSOs were most abundant. These indicators testify that we 
are seeing substantial redistribution of gas, as witnessed by the formation of ionized metal-rich 
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Figure 1: Cosmological density of neutral gas, incidence of CIV absorption, and comoving density 
of luminous QSOs as a function of redshift from QSO absorption-line statistics. Top panel. Mean 
cosmological density of neutral gas, fl 9 , normalized to the critical density (Storrie-Lombardi et al 
1996; Ran et al 1995 (z = 0); Storrie-Lombardi & Wolfe priv. comm.) Bottom panel. Number 
of CIV metal-line absorption systems per unit redshift, n(z) (Steidel 1990); z = 0.3 point from 
Bahcall et al 1993). Filled points from Steidel indicate rest frame equivalent widths IF res t(A1548) > 
0.15 A; open points are for W rest (A1548) > 0.3 A. Hatched areas indicate the range (0 < q Q < 1/2) 
for unevolving cross sections since z = 1.5, beyond which redshift CIV can be measured with 
ground-based telescopes. Bottom panel. Comoving density of optically selected QSOs: filled 
squares from Schmidt et al 1994; open circles from Hewitt et al 1993). H 0 = 50km s _1 Mpc“ 1 , 
qo = 1/2 
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galaxy halos and the efficient fueling of active galactic nuclei. The surge in neutral gas content 
indicates that protogalactic gravitational potentials were deep enough that gas was confined to 
sufficiently high density that it was at least momentarily immune to ionization by star bursts and 
ionizing background radiation. This may be the epoch that disk galaxies formed as secondary 
infall of gas occured into galactic potentials formed in the first round of galactic bulge formation. 
The disk formation would be accompanied by halo enrichment, either by in situ star formation or 
by metal pollution of the extended halo region by winds from the new star forming regions of the 
disk. An alternative view is that this is also likely to be an epoch when small protogalactic lumps 
are merging vigorously, and star burst within the lumps would be effective at ejecting metals into 
an extended region that would at later times constitute the “halo” region of the merger product. 

Figure 1 summarizes only the neutral atomic and ionized gas components. A complete balance 
requires an accounting for all the universe’s baryons, as gas is exchanged between neutral, ionized, 
and molecular phases, as well as the path of stellar evolution leading to the current state where far 
more baryons are contained in stars than in neutral gas. Although the present HI content is only 
~10% of the mass in stars, there was a period at z ~ 2.5 when the HI content apparently was of 
order half of the present stellar mass. Estimates of the mass content in ionized halos suggest that 
they probably contain about ten times the mass contained in the damped Lyman-a absorbers 
at their peak (Petitjean et al 1993). This interpretation of the CIV data relies on theoretical 
modeling, with large uncertainties clue to ionization level and carbon abundance. It is striking 
that the recent HST observations (Bahcall et al 1993) are consistent with no evolution in the 
absorption cross section presented by high column density CIV systems from z ~ 1.3 (Steidel 
1990) to z ~ 0.3, implying that large quantities of ionized gas may still be present, either in the 
form of extended halos or in intergalactic clouds whose mass could far exceed the visible stellar 
mass in galaxies. The hydrogen neutral fraction of these ionized clouds would provide column 
densities of H° well below the regime ordinarily probed by 21cm line observations. At present, 
the integral molecular gas mass content of galaxiers appears to be roughly equal to the neutral 
atomic mass (cf. Kcnnicutt et al 1994). 

Surveys of nearby volumes in the 21 cm line are important in anchoring the z — 0 point 
of Figure 1. A complete inventory of neutral gas in the nearby Universe, of the sort that is 
being provided by unbiased 21cm line surveys, will have tight error bars and thus will carry large 
statistical weight in models that describe the evolution of Q g . Note that the statistical errors 
of the z > 1.6 measurements of are so large that these high z points are consistent with no 
evolution at all. In principle, the low z regime (0 < z < 1.6) of the diagram can be measured 
using QSO absorption line methods in much the same way as the high z points. However, the 
observations are difficult since the Lyman-a line is not redshifted into the optical window until 
z > 1.6, so the low z absorption line work must be done from space observatories such as IUE 
and HST (Lanzetta et al 1995, Rao et al 1995). 

When the data points are plotted as a function of time, as in Figure 2, it is clear that the 
single low z QSO absorption line point applies to a time span well longer than half the age of 
the Universe. Several selection effects make it difficult to obtain a reliable sampling of damped 
Lyman-a absorbers at low z. Cosmological factors, as well as the apparent shrinking of damped 
Lyman-a absorption cross section, with increasing age of the Universe (Lanzetta et al 1995), 
act to make them very rare at recent times: dN/dt oc a{t)t~ k , where k = 2 for q Q = 1/2 and 
k — 3 for q Q = 0. There are also selection effects that are likely to influence the QSO absorption 
line measurements and that affect the low z measurement most strongly: (1) The presence of 
dust in disk galaxies is expected to become increasingly important as they evolve and may act to 
selectively attenuate the light from QSOs behind them, causing these lines of sight to be under 
represented in QSO samples (Fall & Pei 1989, Pei & Fall 1995, Webster et al 1995) although 
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Figure 2: Density of neutral gas as a function of time compared with = 0). £l g (z) from 

references in Figure 1; D stars (z = 0) from Lanzetta et al 1995. Rising dashed curve is KTC model 
for increasing stellar mass; dot-dash is KTC model for declining cold gas content, adjusted as 
described in text to indicate only the atomic fraction. The cross hatched band indicates the range 
of models proposed by Pei and Fall for the true f2 s (z), corrected for selection effects caused by 
dusty damped Lyman-a absorbers. 
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this view is contested (Boyle & Di Matteo 1995). (2) Gravitational lensing may act to selectively 
amplify background objects into QSO samples, but also may bend the light path so as to dodge 
the high HI column densities of the disk (Bartelmann & Loeb 1996, Smette et al 1995a,1996). 

Figure 2 includes curves to indicate trends in stellar evolution relative to the decline in il g with 
time. Recent analyses of the depletion of the integral mass content of galaxies over cosmological 
times due to star formation have been presented by Lanzetta et al (1995) and Pei and Fall (1996). 
A related study by Kennicutt et al (1994, KTC) addresses the prolonging of the current star 
formation rate in z ~ 0 disks due to delayed gas return as stellar populations age. An example 
of the KTC models is presented in Figure 2 to illustrate both the rise of stellar mass with time 
and decline of neutral gas content with time for a disk system without added inflow or allowing 
mass to escape. For this display, the KTC model has been scaled so that the final stellar mass 
and the final HI mass are consistent with the observations at z ~ 0; the relative proportion of HI 
to H 2 has been adjusted for this display to vary linearly with time from Phi/{phi + Ph 2 ) = 1 at 
high z to 2/3 at the present time. The slope of the KTC model at the time marked “Now” in 
Figure 2 can be compared with the steeper slope drawn to indicate the rate at which the current 
star formation rate would consume the present atomic hydrogen content, exhausting the supply in 
only ~1.5 Gyr if there were no additional reservoirs of molecular gas or contribution from delayed 
stellar return. 

The models of Lanzetta et al require a balance between the decline of Q g and a rise in stellar 
mass; along with assumptions that stellar return occurs instantaneously as each generation of stars 
is formed, this constrains the mean star formation rate history of the Universe. A consequence 
is that their models result in an uncomfortably large number of stars at z ~ 0 with low metal 
abundances. Pei and Fall suggest that this problem can be solved with a family of dusty models, 
which imply that the damped Lyman-a statistics drastically underestimate the integral gas content 
at all redshifts. The current generation of radio telescopes are not sensitive enough to resolve this 
controversy by simply looking back to measure the integral HI content at z ~ 1/2. On the other 
hand, choosing complete samples of radio selected high z quasars for background probes would 
remove possible selection effects by dust. 

At z = 0, recent 21 cm line measurements are indicating that the bulk of the atomic hydrogen 
content of the nearby universe is bound into galaxies with optical counterparts (Zwaan et al 1996, 
Schneider 1996, Szomoru et al 1994, Henning 1995, Briggs 1990). Furthermore the normalization 
of seems to be well understood (cf. Zwaan et al 1996, Rao & Briggs 1993), although there is 
still concern over the normalization of even the optically determined luminosity function (cf Ellis 
et al 1996, Glazebrook et al 1995). Clearly the determination of the integral HI content is a 
measurement that the Parkes Multibeam Survey will clarify since it will be complete, unbiased by 
extinction and optical surface brightness, and will have well understood sensitivity limits. 

2 Distribution of N(H°) Column Densities 

Another place where future surveys can play an important role will be in exploring lower HI column 
densities than have been observed in the 21 cm line in the past. QSO absorption line statistics 
over a wide range of Lyman-a line strengths specify that the incidence of absorption becomes 
increasingly prevalent toward lower column densities, so that along a randomly chosen sight line, 
the probability of interception rises by roughly a factor of 10 for every decrease by a factor of 100 
in Nhi■ This behavior is quantified by the f(Nni) distribution, defined so that f(NHi)dXdNHi is 
the number of lines detected within a range dNm centered on N HJ over a “normalized absorption 
distance” dX, where X = |[(1 + z ) 2 — 1] is the normalized absorption distance from zero redshift 
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to z for q 0 = 0. The function f(N Hr ) is roughly proportional to X -1,5 over the range Nhi = 10 13 
to 10 22 cm -2 , although there is evidence for subtle structure possibly related to opacity in the 
Lyman continuum that occurs for layers with Nhi > 3xl0 17 cm -2 (Petitjean et al 1993). At 
high redshift, the frequency of absorption for these optically thick absorbers is ~0.9 per unit X 
(Petitjean et al 1993); for Nhi > 10 18 cm -2 , the frequency is roughly halved. The derivation of the 
/( Nhi) at these Nhi is especially uncertain, since the entire Lyman series is heavily saturated in 
the regime where the Lyman continuum is optically thick, and the damping wings, which permit 
an unambiguous measure of column density, do not become readily observable until Nhi is well 
in excess of 10 20 cm -2 . Thus, column density measurements in this regime are very uncertain, 
leading Petitjean et al to plot only one point on their f(Nni) diagram for 10 177 to 10 20 ' 5 cm -2 . 

Recent large HI surveys with filled aperture telescopes are routinely capable of detecting column 
densities below 10 19 cm -2 (Schneider 1996, Zwaan et al 1996, Briggs et al 1996), provided the 
emission fills the telescope beam. The Arecibo survey by Sorar (1994; see also Briggs & Sorar 
1996) was optimized to be sensitive to Nhi = 10 18 cm -2 (5<r), and the survey observed over 5000 
independent beam areas to a depth of 7500 km s -1 , covering a total absorption path AX « 120. 
The Arecibo beam subtends ~3 Kpc at 3 Mpc and 70 Kpc at 75 Mpc, which is a reasonable match 
to the cloud sizes deduced for the Lyman-a forest (cf. Smette et al 1995b, 1992). To date, only 
one of the 61 detections (Zwaan et al 1996) has not been identified with a high column density 
layer of the type associated with the neutral intergalactic medium of a galaxy, implying that a 
separate population of low column density objects can add only a small fraction of the current 
integral HI content already identified with galaxies. 

The high redshift f(N H i) distribution would imply of order 50 interceptions in the range 10 18 
to 10 19 cm -2 for the pathlength X explored by the Arecibo survey. Where are they? At least a part 
of the discrepancy is likely clue to evolution of Lyman-a forest cloud population. An additional 
observational problem is that the high column density end of the forest cloud distribution (around 
10 1 ' cm -2 and above) has associated metal lines, such as CIV and Mgll, which has historically 
caused then to be identified with hypothetical galaxy halos; single dish observations seldom have 
the resolution to reliably separate the HI signal from a halo of a spiral galaxy from the bright 
signal originating in the main body of the galaxy, unless there are strong kinematic effects that 
create a difference in gas velocity as a function of radius and the halo gas is very extended. 

Further considerations in the study of this intermediate column density range are the theoret¬ 
ical models that consider ionization of extended gas around galaxies by the extragalactic ionizing 
background (Sunyaev 1969, Corbclli et al 1989, Maloney 1993, Charlton et al 1994). Many of these 
models predict a strong dip in f(Nni) between 10 17 ' 5 and 10 19 ' 5 cm -2 . If similar arguments apply 
to a population of intergalactic clouds or super-LSB galaxies, then an interesting experiment now 
coming into the realm of possibility will be to push the sensitivity limits of the local HI emission 
observations down to ~10 1 ‘ cm -2 where ionized layers of high column density might be detected. 

3 Where to find “HI Clouds” Now 

Some very interesting examples of HI without coincident optical emission have turned up by chance 
in “off-scans” from 21cm line studies (Schneider 1989, Giovanelli & Haynes 1989, Chengalur et al 
1995, Giovanelli et al 1995), although all of these detections appear to be associated in some way 
with nearby visible galaxies. For example, the protogalaxy of Giovanelli & Haynes (1989) now 
appears in VLA observations to more closely resemble a galaxy with tidal remnant (Chengalur 
et al 1995) than a single large cloud. An interesting consequence of high sensitivity VLA surveys 
of HII galaxies and LSB dwarfs (Taylor et al 1993; Taylor et al 1996) is an apparent increased 
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probability of finding dim Hi-rich companions to the HII galaxies as compared with the LSBs. 
This kind of study will not be possible with the crude resolution of a single-dish surveys. The 
tendency of small galaxies to be found in the vicinity of large ones rather than in isolated regions 
(cf. Szomoru et al 1996) will complicate the derivation of a Hi-mass function that extends to faint 
masses when survey observations are made with a large beam. On the other hand, isolated HI 
clouds (if they exist) and HI clouds associated with early type galaxies, such as polar rings and 
fresh examples of the Leo Ring (Schneider 1989), should be readily identified if confusion with 
other galaxies in the same field is not too much of a problem. 

4 Conclusion 

Large sky surveys in the 21cm line will improve the measurement of the integral HI content of 
the Universe, and this value will find immediate use in anchoring theories of galaxy evolution. 
Telescopes with large beams will suffer from confusion, which will complicate the determination 
of the faint end slope of the Hi-mass function. The low column density regime with N HJ < 10 18 
cm -2 is a new frontier that awaits exploration in the local universe. 
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